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Abstract 


The icy Galilean satellites of Jupiter — Europa, Ganymede, and Callisto — have unusual radar scat- 
tering properties compared with those of the terrestrial planets or Earth’s Moon. There are three 
main features of the data that distinguish these targets: (1) The radar cross-section normalized by 
the geometrical cross-section is an order of magnitude larger than that of any terrestrial planet. (2) 
The reflected power is almost evenly distributed between two orthogonal polarizations with more 
power being returned in the same circular polarization as was transmitted whereas virtually all of 
the power returned from the terrestrial planets is contained in the opposite circular polarization to 
the one that was transmitted. (3) The echo power spectra have a broad shape indicating a nearly 
uniformly radar-bright surface in contrast to the spectra from the terrestrial planets that contain a 
strong quasi-specular component from the vicinity of the sub-radar point and very little reflected 
power from the rest of the surface. The normalized radar cross-sections decrease as the areal water- 
ice coverage decreases from Europa to Ganymede to Callisto. Recently, radar echoes from the polar 
caps of Mars and Mercury, and from Saturn’s satellite Titan imply similarly strong cross-sections 
and have classically unexpected polarization properties and it is also thought that this is due to the 
presence of ice on the surface. 

We analyze a model called the radar glory model and show that the main features of the radar 
echoes calculated from this model agree well with the observations from all three icy Galilean satel- 
lites. This model involves long radar paths in the ice below the surface and special structures in 
which the refractive index decreases abruptly at a hemispherical boundary. It is not known whether 
such structures exist or how they could be created, but possible scenarios can be imagined such as 
the formation of an impact crater followed by deposition of a frost layer follwed by a resurfacing 
event in which a layer of solid ice is placed above the layer of frost. Regardless of the exact geophys- 
ical processes required to create such structures, the superior ability of this model to account for all 
of the important observations with very few adjustable parameters and with no ad hoc assumptions 
is a compelling argument in support of at least the electromagnetic model. The key features of the 
electromagnetic model are multiple subsurface scattering events, total internal reflection, and a low 
degree of randomness imposed on a deterministic geometry that strongly favors backscattering. 
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Chapter 1 


Introduction 


The field of radar astronomy began shortly after World War II when radar echoes were obtained 
from Moon. (For a brief history of radar astronomy see Ostro (1987),Muhleman et al. (1965), and 
Muhleman et al. (1965).) In the following few decades before the echoes were obtained from 
the Galilean satellites, echoes were obtained from Mercury, Venus, Mars, Saturn’s rings, and a few 
asteroids. The echoes from the terrestrial planets and Moon are similar and set the standard for what 
was expected from planetary targets. The typical echo from the terrestrial planets is characterized by 
a small radar cross-section relative to its geometrical cross section. When circularly polarized fields 
of a given helicity are transmitted nearly all of the reflected fields are in the opposite helicity state 
and the spectra (a plot of power versus Doppler frequency) show a strong quasi-specular component 
around the sub-radar point and a weak diffuse component from the rest of the surface. 

The first attempt to detect radar echoes from the Galilean satellites was made in 1970 (Campbell 
et a/., 1977), but was not successful. In 1974 Ganymede became the first Galilean satellite detected 
by radar (Goldstein and Morris, 1975) using the 64 meter Goldstone antenna. From the very weak 
echoes obtained it was estimated that the radar cross section of Ganymede was 12% that of a 
perfectly conducting, smooth sphere — a result that was comparable with other radar targets such 
as the terrestrial planets and asteroids. These echoes also revealed a broad power spectrum which 
was quite different than that observed for the terrestrial planets, but comparable to the spectrum 
obtained from an asteroid’s rough surface. This detection was of interest for the fact that it could 
be done, but did not hint at the truly unique nature of Ganymede’s, as well as the other icy Galilean 
satellites’, radar echoes. The entire set of Galilean satellites were observed using radar in 1975 
using the Arecibo S-band radar system which had five times the sensitivity of any other system for 
studying the Galilean satellites (Campbell et a/., 1977). Since then, radar echoes from the Galilean 
satellites have been obtained whenever Jupiter was at declinations favorable for viewing by Arecibo 
(a period of about four years once every about eleven years). A summary of the observations made 
in the years 1975-1979 are given in Ostro (1982) and a summary of the observations made in the 
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years 1987-1991 are given in Ostro et al. (1992). 

The radar echoes from the Galilean satellites are very different than those from the terrestrial 
planets and the term bizarre has often been used when referring to them. Soon after it was recognized 
that the echoes from these satellites were so unusual it was realized that the fact that the surfaces are 
made up of large quantities of water ice, which is a poor absorber of microwaves, had something to 
do with the observations. However, after more than almost two decades in which at least five serious 
models for explaining the data have been proposed the issue has not been settled to any degree of 
certainty. The basic reason for this is that the data do not sufficiently constrain the models. New 
types of data are required such as data from a bistatic radar experiment in which different antennas 
at different angular positions are used to transmit and receive the radio waves so that the angular 
dependence of the scattering can be determined. 

In the following chapter we will give background information on the basic physical characteristics 
of these satellites, the properties of ice that are important to understanding the radar echoes, the 
experiment and the data, and the models that have been proposed to explain the data. Then in 
chapter 3 we present the details of our analysis of the radar glory model. This model was first 
proposed by Eshleman (1986) and this report extends the analysis of that model to include certain 
effects that were previously neglected and to calculate the Doppler spectra which had not been 
calculated before. That chapter will end with a comparison of model calculations with the data and 
it will be seen that the model predictions are consistent with all the main features in the data. 



Chapter 2 


Background 


2.1 The Incidence and Importance of Ice 

Four of the satellites of Jupiter are referred to as the Galilean satellites in honor of Galileo who 
is credited with their discovery in 1610 using one of the earliest telescopes. Of these four three 
of them — Europa, Ganymede, and Callisto— are different than the fourth one — Io — because their 
surfaces are largely composed of ice (throughout this document the term ice refers to H 2 O ice, 
usually in its main crystalline form called ice Ih (Hobbs, 1974)). It is this distinguishing feature that 
is most likely responsible for the unusual radar backscattering properties of the three icy satellites 
that is the subject of this report. Io is a very interesting satellite, but it is not a very good radar 
backscatterer and so I will not discuss Io any further. 

Some of the physical characteristics of the Icy Galilean satellites are summarized in Table 1. 
Europa is about the size of Earth’s Moon and its bulk composition is mostly silicate rock material 
with a crust of ice to a depth of at least 25 km. Spectral measurements indicate that the surface of 
Europa is virtually 100% ice. Ganymede and Callisto are about the size of Mercury and their bulk 
composition is a mixture of ice and rock. Ganymede is a well differentiatied body and so much of 
the ice is contained in the upper layers. Spectral measurements indicate that Ganymede’s surface 
is about 50% ice. Callisto on the other hand is not very well differentiated and so it does not have 
most of its ice near its surface. Spectral measurements indicate that its surface contains less ice 
than Ganymede and that the aerial coverage is about 10%. The temperatures of all these satellites 
is about 100 K. 

Most of the theories that have been proposed to explain the strange backscattering properties of 
these satellites involves subsurface scattering and long radar paths. This is possible because ice at 
these cold temperatures is not a strong absorber at microwave wavelengths. Actually no laboratory 
measurements are available at these temperatures of the complex dielectric constant. It is usually 
assumed that the absorption length in ice on these satellites is much longer than what is required for 
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the particular model being analyzed. This assumption appears to be based on an extrapolation using 
the standard Debye theory of the permitivity with the temperature dependence of certain parameters 
determined from theory and measurements made at long wavelengths (A > 1.5 km) and temperatures 
down to 200 K (Thompson and Squyres, 1990). However, Thompson and Squyres (1990) argue 
that this extrapolation is inaccurate at microwave frequencies because the Debye formulas don’t 
account for the large tails of absorption features in the millimeter wavelength range that influence 
the microwave absorption. Thompson and Squyres derive formulas with empirically determined 
parameters based on measurements from optical wavelengths to 8.6 m-A at temperatures above 213 K 
which they propose to be better than the previous formulas for extrapolating to cold temperatures 
at microwave frequencies. The real part of the dielectric constant given by the Thompson and 
Squyres formulas agrees with the previous results at all wavelengths and temperatures. However, 
the imaginary part of the dielectric constant in the Thompson and Squyres formulas is orders of 
magnitude larger than that previously reported. A convenient parameter for discussing absorption is 
the absorption coefficient which is proportional to the imaginary part of the dielectric constant and 
is the propagation distance required to reduce the power by a factor 1/e. At 12.6 cm-A and 100 K, 
the Thompson and Squyres formulas give an absorption length of 22 km whereas the older formulas 
give 9 x 10 15 km. The absorption lengths at all temperatures of interest at the Galilean satellites 
based on the standard formulas are much larger than the radii of the satellites. This justifies the 
usual analyses that regard the ice to be free of absorption. However, if the Thompson and Squyres 
formulas are correct, then the absorption in the ice may play an important role in the differences 
in the observed scattering properties of these satellites. For future reference, the refractive index 
and absorption length based on the Thompson and Squyres formulas at temperatures of interest 
are tabulated in Table 2. In our analysis, we have neglected the effects of absorption by the ice 
like all the other theories. The effects of absorption can be accounted for in an ad hoc manner by 
multiplying the power by exp( — L/c*) where L is the total distance traveled in ice by the radar waves 
and a is the absorption length. 

Table 3 shows the possible polymorphs of ice that could possibly exist on the icy Galilean satellites 
(Gaffney and Matson, 1980). The polymorphs are forms of ice other than the normal Ih hexagonal 
form of ice that form at higher pressures. Gaffney and Matson hypothesize that the polymorphs 
will occur when impacts occur during crater formation. They also expect that at the temperatures 
found in the outer solar system, that the I c form is stable and will occur in abundance since there 
has been a long time for it to accumulate. A likely place to find this form, they claim, might be as 
a light frost cover in the low albedo regions of Ganymede poleward of 40° latitude. 

Meteor bombardment, pressure, and thermal effects also act to modify the structure of ices on 
satellites (Smoluchowski and McWilliam, 1984). Meteor bombardment in porous ice results in some 
densification, but mainly results in the formation of ice crystalline polymorphs and in amorphous 
ice below 150 K. Smoluchowski (1983) discusses amorphous ice which is created when water vapor 
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condenses directly into ice at temperatures below about 150 K. Smoluchowski hypothesizes that 
there should be up to 50% of the ice layer at an impact site which is amorphous ice. Impacts of 
meteorites into ice may produce crystalline ice polymorphs in the crater; the ejecta will consist of 
both liquid and vapor; the liquid will partially evaporate and the rest will refreeze on the surface 
forming regular crystalline ice I h ; the vapor will partially escape from the planet and partly condense 
into amorphous ice on the solid ejecta and around the crater. Pressure densifi cation, controlled by 
creep or slow plastic deformation and diffusion, is significant to depths of hundreds of kilometers 
because of the low temperatures involved. Thermal effects result in densification as well due to: 
1) migration of pores toward warmer regions when thermal gradients exist; and 2) diffusion when 
isothermal conditions exist. 

There doesn't appear to be much data on the dielectric properties of these other forms of ice, 
particularly at the cold temperatures of interest for the Galilean satellites, but it is likely that the 
dielectric constant is differnet and this may cause interesting inhomogeneities and discontinuities in 
the refractive index of the ice. 


2.2 The Unusual Radar Backscattering Properties 

Figure 1 schematically shows the geometry of a radar backscattering observation of the Galilean 
satellites using Arecibo. The optimal viewing times are at Jupiter opposition (when Jupiter, Earth 
and Sun are in alignment with Jupiter at its closest to Earth) when Jupiter is about 4 AU (astro- 
nomical units) from Earth. The round trip light time is a little over an hour. In each experimental 
“run” CW (continuous wave, meaning a sinusoidal wave that is unmodulated for a time very long 
compared to the period of oscillation) signals were transmitted for the duration of the round trip 
light time when the radar was switched to receive mode. Usually circularly polarized fields in one 
helicity state were transmitted and the fields in both helicity states were received. On a few occa- 
sions linear polarization in one particular orientation relative to the feed of the Arecibo antenna at 
transmit time were used and two orthogonal linear polarizations were received. 

The use of CW signals permits a certain kind of mapping of the surface of a satellite as described 
in Green (1968) (see also Ostro (1987)). Figure 2 shows a satellite of radius R and spin frequency Q 
in radians. A cartesian coordinate system with its Z axis pointing in the direction of Earth and with 
the XZ-plane containing the spin axis of the satellite is defined. Also shown is a scattering center 
located at spherical polar coordinates (i£, 0,<$). If the incident fields are sinusoids of frequency 
u>o, then the fields reflected from the scattering center will have a doppler shifted frequency. The 
Doppler shifts incurred in the signal received by the radar at Earth are due to many sources such as 
the orbital and rotational motions of the Earth, the orbital and rotational motions of the satellite, 
and other less important contributions such as Gravitational red and blue shifts. The radar system 
at Arecibo is programmed to remove the Doppler shifts due to Earth’s motions and the orbital 
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motions of the satellite; these motions are known very accurately. The remaining Doppler shift is 
due to the spin of the satellite. If we call u the frequency of the reflected fields after these known 
effects are removed, then it is a simple matter to calculate 

w — wo = 2koY s Q, sin ^ 

where Jbo = u>o /c is the wavenumber of the incident fields and Y s is the Y coordinate of the scattering 
center. From this formula we see that if we compute the Fourier spectrum of the signal that is 
received at Earth, then the power at frequency u> has contributions from all scatterers of the same 
Y coordinate across the surface of the satellite. Thus, if we look at the power at a single frequency 
it is as if we looked at the power reflected in a strip parallel to the plane containing the spin axis 
and the radar line of sight. 

Another possible mode of operation of the radar that would provide a different sort of map would 
be to transmit short pulses. Then the power reflected at a given time before the arrival of another 
pulse is due to scattering centers located at the same colatitude 0. Thus, this would provide a kind 
of map of the surface in terms of rings of constant colatitude. It is also possible to transmit sinusoidal 
signals that are modulated by short pulses which would then allow the power to be resolved in time 
and frequency. The power at a specified time and frequency would be from scattering centers at the 
intersection of a ring of constant 0 and a line of constant V, that is, from two points on the surface. 
It turns out that it isn’t possible to do this kind of mapping on the Galilean satellites as described 
here. The reason is that in order to prevent aliasing of the Fourier spectrum, the frequency of the 
pulses would have to be larger than the bandwidth of the reflected signals and in order to prevent 
a kind of aliasing in the time dispersed signal (to arrange it so that only one pulse is interacting 
with the surface of the satellite at a time) it is necessary that the time between the pulses be larger 
than the round trip travel time of a pulse from the sub-radar point (at 0 = 0) to the terminator 
(at 0 = 7t/2) and back. Thus, these two requirements place constraints on the time between the 
pulses, but the constraints are in opposite directions. If B is the bandwidth of the reflected signals 
and if T is the round trip time of a pulse across the surface and if r is the time between pulses, then 
the two requirements are 1/r > B and r > T which can be combined to one requirement for the 
combined type of mapping 1 > BT. Assuming = 7r/2, then B — 8tt R/(\qP) where Ao = u>o fc 
is the wavelength of the incident fields and P = 27rR/Q is the rotation period of the satellite and 
T = 2 R/c. Table 3 tabulates B , T } and their product. From the table we see that the product 
BT is larger than unity in all cases except Callisto at 70 cm-A which is not a good wavelength in 
terms of SNR; therefore, it is not possible to simultaneously map the surface in frequency and time 
delay. Radar targets for which BT is greater than unity are said to be overspread. All of the radar 
echoes from the Galilean satellites have been obtained using CW signals. Recently Hudson and 
Ostro (1990) developed a method to produce two dimensional maps of the surfaces of spherically 
symmetric targets using only Doppler shifted echoes from CW signals from several different rotation 
phases (that is, from several different runs of the experiment on different nights). The results of this 



CHAPTER 2. BACKGROUND 


7 


mapping procedure using the most recent Arecibo data are discussed in Ostro et al. (1992). 

The analysis presented here deals with the basic mechanisms causing the unusual backscattering 
properties of these satellites. Therefore, I will not discuss mapping techniques or correlations between 
radar results and optical images. 

Figure 3 shows a comparison between a typical radar echo from a terrestrial planet and from 
Europa. The main features of this spectrum from Europa is typical of what is also observed from 
Ganymede and Callisto and will serve for the purpose of this discussion as a typical spectrum from 
the Galilean satellites. Other examples of spectra from Ganymede and Callisto will be shown later. 
Table 5 also summarizes the main features of all of the radar echoes obtained up to 1982 (Ostro, 
1982). For comparison, numbers are also given for Earth’s Moon which is typical of those for the 
terrestrial planets; Moon has been studied the most using radar. The first thing that is obviously 
different between the icy Galilean satellites spectra and that of a terrestrial planet is the size of 
the spectra. The spectra shown in Figure 3 are normalized by the geometric cross-section of the 
respective radar target. The total area under the spectrum from the icy Galilean satellite is much 
larger than that of the terrestrial planet as can also be seen from the first row of Table 5. Europa 
and Ganymede have larger radar cross-sections than would a perfectly conducting sphere of the same 
radius as these bodies which would have <7=1. Another striking difference is that in the Europa 
spectrum, the spectrum marked SC, which is the spectrum observed in the same circular polarization 
helicity state as was transmitted, is larger than the spectrum marked OC, which is the orthogonad 
circular polarization helicity state. This is the reverse of the situation that is observed from a 
terrestrial planet. Although the linear polarization ratio is not larger than unity, it is significantly 
larger than that observed for terrestrial planets. The final difference that is of significance for the 
work here is that the spectrum from the terrestrial planet has a strong feature near zero Doppler 
frequency, the quasi-specular component, and small refelctivity at other Doppler frequencies, the so 
called diffuse component. The spectra from the icy Galilean satellites don’t appear to divide into 
these two components; they have significant power at all Doppler frequencies and they don’t have 
a strong, sharply peaked feature at zero Doppler frequency. The power law exponent, m, shown in 
Table 5 for the Galilean satellites indicates a broad power spectrum as in Figure 3. If every point 
on the surface scattered isotropically over a hemisphere, then the exponent would be 2; this case is 
referred to as a Lambert scattering law. The exponent for Moon refers only to the quasi-specular 
component. The exponent for the diffuse component is probably closer to the size of that of the 
Galilean satellites. 

Recent radar investigations of Saturn’s satellite Titan (Muhleman et ai } 1990) and the polar 
ice cap of Mars (Muhleman et al 1991) show similar radar properties to those of the icy Galilean 
satellites. 
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2.3 Previous Models and Analysis 

In the next chapter I will discuss the analysis I have done of a model first proposed by Eshleman 
(Eshleman, 1986a; Eshleman, 1986b). In this section I will discuss four other models that have been 
proposed to explain the unusual backscattering properties as well as the work done by Eshleman 
prior to the work to be discussed in the next chapter. Figure 4 shows the essential geophysical 
structure involved in the five models. Three of these (Ostro and Pettengill, 1978; Goldstein and 
Green, 1980; Hagfors et a/., 1985) preceded the radar glory model of Eshleman and this model 
shares several of the features of the others. One model (Hapke, 1990) appeared after Eshleman ’s 
model. It is perhaps the most attractive model in terms of its geophysical plausibility, but there is 
no way to use the data in hand to discriminate between the different models. Other types of data 
such as bistatic radar experiments to explore the phase function of the scattering centers are needed 
to completely understand the mechanisms causing the unusual scattering by these satellites. 

Most of the models involve multiple scattering in order to reproduce the observed polarization 
behavior. However, the radar glory model of Eshleman and the surface crater model of Ostro and 
Pettengill involve a small degree of randomness in the placement of the scatterers in a basically 
deterministic and special geometry whereas the models of Goldstein and Green and Hapke involve 
complete randomness in the location and orientation of scatterers. The model of Hagfors et ai is 
the only one that does not involve multiple scattering and doesn’t invoke any randomness in its 
inital conception. Eshleman (1986), however, showed that this model would predict a null in the 
scattering in the backscattering direction unless there were some degree of incoherence involved. 

Ostro and Pettengill (1978) showed that the observed polarization properties were consistent 
with a model in which the radio waves are reflected twice by hemispherical craters on the surface. 
However, since the craters are on the surface, much of the power is transmitted into the ice and this 
model is not able to account for the huge radar cross-sections without making the unlikely assumption 
that the surface refractive index is closer to that of silicate rocks (n = 2.4) than that of ice (n = 1.8). 
This is especially unlikely since it would require Europa to have its surface refractive index be the 
closest to that of silicate rocks and it is known that Europa’s surface, of the three icy Galilean 
satellites, contains the highest ice to rock ratio (see Section 2.1). Ostro (1982) concluded that this 
model was not very plausible because of the high refractive index material required and because it 
requires craters that are perfectly hemispherical and too smooth on scales of the wavelength. 

Goldstein and Green (1980) proposed that the upper few meters of the surfaces of these satellites 
are full of randomly oriented cracks and fissures that form voids in the ice with planar facet interfaces 
between the ice and vaccuum. They recognized that when a plane wave is incident on a planar facet 
at angles greater than the Brewster angle that the incident sense of polarization is largely preserved. 
They ran Monte Carlo simulations of photons reflecting from a slab of ice containing such voids 
with the orientation of the normal to the planar interfaces being completely random and found that 
this model could account for the observed poarization properties of the icy Galilean satellite radar 
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echoes. In the best-fit conditions, most of the photons were totally internally reflected from the 
voids. They found that the scattered energy behaved very much as if it had been scattered by an 
ideal diffusing screen; thus, the direction of the energy is isotropic. Ostro (1982) found that this 
calculation of Goldstein and Green had a slight error in it that overestimated the ability of the model 
to invert the usual polarization behavior and that the model could not simultaneously be fixed to 
fit the polarization behavior and the huge radar cross-sections. Ostro argues that other scattering 
mechanisms such as scattering by wavelength sized particles could help to increase the radar cross- 
sections, but that these single scattering events would rapidly destroy the ability of the model to 
invert the polarization. Nevertheless, he is able to find a match to the data by requiring that the 
ice refractive index be very small (n = 1.2 in the case of Europa). He points out that such low 
refractive indices imply low mechanical strength and hence make the assumption that voids could 
be maintained untenable. He salvages the model by proposing that the voids could be replaced by 
material of slightly different refractive index from that of ice and of similar strength. 

Hagfors et ai (1985) concluding that the Goldstein and Green model can’t account for the huge 
radar cross-sections propose a model involving refraction by spherically symmetric inhomogeneities 
in refractive index below the surface of the ice. They assume that the inhomogeneities are large 
enough that geometrical optics is valid. Since there are no reflections involved and the incident 
fields are simply turned around toward the radar, the incident polarization is preserved entirely. 
Also, they find that there is a caustic in the backscattering direction and thus in the geometrical 
optics approximation the field is infinite. Thus, they conclude, without calculating the actual radar 
cross-sections implied by their model which can’t be calculated using geometrical optics, that if this 
scattering mechanism operates in conjunction with other mechanisms that the radar properties of 
the Galilean satellites are explained by it. 

Eshleman (1986) computes the wave optical correction to the infinities predicted by geometrical 
optics in the Hagfors et ai model and finds that when polarization is taken into account that the 
backscattering from perfectly spherical backscattering lenses would be zero. The principal contribu- 
tion to the backscattered field comes from a ring or halo around the scattering center and is similar 
to the optical effect called the glory in which such halos are observed around the shadow of an 
observer when the incident light comes from behind the observer and the shadow is cast on water 
drops in a cloud. To salvage this model, he was lead to propose that the glory halos broke up into 
patches or glints due to azimuthal imperfections in the shape of the lense. The fields were coherent 
within each such glint that were assumed to be larger than a wavelength in size and the fields from 
different glints added incoherently. 

Eshleman (1986) proposes a modification to the model in which the inhomogeneity in refractive 
index is a step discontinuity at the wall of a buried crater in which the refractive index goes from 
some value above the wall to a smaller value below. He argues that there are a number of ways in 
which such structures could occur. One way involves a cratering event followed by a process (either 
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sublimation or frost deposition) whereby the surface ice becomes less dense and hence its refractive 
index becomes smaller than the more solid ice below which is then followed by a resurfacing in which 
the frost layer is overlain by a more dense ice overburden. Another possibility is that the cratering 
event itself compacts the ice near the surface of the crater walls which results in a layer of higher 
refractive index material lining the walls or perhaps the rocky material of the impacting object lines 
the walls of the crater. In any case, the main features of the model are a hemispheroidal geometry 
that strongly favors backscattering over scattering in any other direction, a decrease in refractive 
index at the hemispherical boundary that permits multiple total internal reflections thereby reflecting 
all the incident power at each reflection and also causing a differential phase insertion in two modes 
(TE and TM) of reflection, and a slight randomization of the shape that breaks up the glory halo 
into several coherent patches. This can be taken purely as a prescription of the electromagnetic 
model without being too specific about the mechanical model that goes with it. The key feature 
distinguishing this model from the Goldstein and Green model and the Hapke model (to be discussed 
next) is the specialized structure and the relatively low order of randomness involved in the Eshleman 
model. The result is that the Eshleman model easily accounts for the polarization inversion and the 
large cross-sections observed. 

Eshleman ’s analysis does not account for the effects introduced by the interface between the 
overburden and the vacuum above and he doesn’t calculate the shape of the Doppler spectra from 
this model. The work reported in Chapter 3 completes the analysis of the radar glory model. In 
that chapter we also consider a possible explanation for the slight frequency dependence indicated 
in Table 5. 

Hapke (1990) and Hapke and Blewett (1991) propose a different model which is like the Goldstein 
and Green model in that it involves mutiple scattering below the ice surface, but it does not involve 
total internal refelctions from the scatterers. Hapke proposes that the scatterers are silicate rocks 
of about wavelength size suspended in the ice. He points out that despite the randomness involved 
there is a coherence involved in the backscattering direction because for any photon which can find 
its way through the multiple scatterings back to the radar there is another one that was transmitted 
from the radar that found its way along the same path but in the reverse direction. Thus those two 
photons will add coherently. The result is that there is an enhancement in a small range of angles 
about the backscattering direction where the scattered intensity rises from a diffuse background 
level to values as much as two times this level in the backscattering direction. He calls this effect 
the coherent backscattering opposition effect (CBOE). This doubling of the power would also occur 
in the Goldstein and Green model as well as the Eshleman model. It doesn’t appear to have been 
included in the Goldstein and Green model, but it has been included in the Eshleman model although 
he didn’t point it out explicitly presumably because only backscattering results were being discussed 
and the only significance of the factor of two for this one direction is that the solution of the problem 
is incorrect without it. (See section Chapter 3 where the factor of two is explicitly pointed out and 
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the resulting formulas can be easily seen to agree with Eshleman’s.) Hapke does a scalar radiative 
transfer analysis and hence his calculation of polarization properties is done in an ad hoc manner. 
Nevertheless, Hapke concludes that his model predicts the observed polarization properties, the huge 
cross-sections, and the broad Doppler spectra. However, to obtain a fit to the large cross-sections he 
requires scatterers with albedos approaching untiy, that is, with low loss. Peters (1992) did a more 
rigorous perturbation calculation including polarization effects for the Hapke physical model and 
found that if the scatterers are anisotropic, that is, they are forward scatterers implying they are 
of wavelength or larger size, then the polarization ratios computed from his analysis are consistent 
with those observed. In his paper he only does the comparisons with the polarization ratios and 
doesn’t discuss conditions necessary to fit the cross-sections. 


Chapter 3 


Analysis of Radar Glory 
Backscattering 


3.1 General Considerations for Scattering by Satellites 

We consider plane waves incident on a spherical planet as shown in Figure 2. We define a satellite- 
centered cartesian coordinate system with unit basis vectors (e x ,e yi e z ) having e z pointing toward 
Earth and e x , e z , coplanar with the spin axis, J?, of the planet. We will consider a scattering 
center located at the point (X,Y,Z). It will be convenient to also define angular coordinates (0,$) 
of the scattering center as shown in Figure 2. 

The incident field is a time harmonic travelling wave with frequency u 0} wavelength A 0 , and 
wavenumber k 0 = 27t/Ao. The medium is vacuum with permeability /i o, and permittivity eo so 
that the wave phase velocity is c = y/^io^o, and the wave impedance is t) q = \J f^o/so. In the 
planet-centered coordinate system the incident field is, 

E i = (E i x e x +E' Y e Y )e ( - ikZ - iu ° t) (1) 

H' = — — ez x E' (2) 

m 

where E l x and E l y are constant, complex amplitudes. We will consider both linearly polarized and 
right-hand circularly polarized incident fields. The linear polarized case will have its electric field 
vector at an angle 4>p with respect to the X-axis. For the linearly polarized case, 

E l x = £*cos<l>p (3) 

E l Y = E i sin <I>p . (4) 

For the circularly polarized case, 

E' x = E i (5) 


12 
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E\ = -if? 1 ’ (6) 

We are interested in calculating the radar cross-section for backscattering of two orthogonal 
components of the electromagnetic field, either the same linear (SL) and orthogonal linear (OL) 
or the same circular (SC) and orthogonal circular (OC) relative to the transmitted polarization. 
The reflected field, Ep , where P indicates the polarization, is time varying with a deterministic 
component that is a rapidly oscillating signal modulated by the Doppler shifted and time delayed 
contributions from the many scattering centers and a random component due to the many scattering 
events that contribute to the backscattering from a particular scattering center which is usually 
treated in a statistical manner. (Actually, there is no randomness involved because there is only 
one realization of the scatterers and there is no relative motion, random or otherwise, between the 
scatterers. Thus, the detailed time variation of this ‘random’ component of the reflected signal 
contains information about the scattering process, but it is too much information relating to the 
properties and locations of each individual scatterer. There is no way to use all this information. 
In order to reduce the amount of information we forgo any detailed knowledge of the properties 
and locations of the individual scatterers and opt for a statistical treatment of those properties in 
which we only describe the mean and perhaps standard deviation of those properties and hence of 
the time varying reflected signal.) Define Ep(X,Y) to be the backscattered field received at Earth 
due to a scattering center of area dA located at the point (X, Y, V ft 2 — X 2 — Y^) on the surface of 
the satellite (refer to Figure 2). This field is time varying consisting of a time harmonic component 
at the Doppler shifted frequency u> given by 


u> — wo = 2fco YQ sin W 


( 7 ) 


and a random component due to the many backscattering events included in the area dA (neglecting 
other time variations that are removed as discussed in Section 2.2). The radar cross-section of this 
scattering center is then defined as (Green, 1968), 


a P {X, Y) - 4nr E ^ 

where the line over \Ep(X, Y)| 2 indicates an infinite time average. I will write this as 


( 8 ) 


ap(X ) Y) = N{a P (X,Y)) (9) 

where (^(AjY)) indicates the expected value (t. e., ensemble average) of the radar cross-section 
for a single backscattering event and N is the number of backscattering events within the area dA. 
We also define the specific cross-section (the cross-section per unit area of surface also known as the 
surface scattering coefficient (Ulaby ei al. , 1981)), 

a° P (X, Y) = ^- = p(X, Y)(<r P (X, Y)) 


( 10 ) 
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where <jp is the total cross-section from all scattering centers and p(X,Y) is the number density of 
backscattering events on the surface. 

The total reflected field Ep is again time varying due to the time delays and Doppler shifts from 
the individual scattering centers. Since the icy Galilean satellites are overspread and CW signals are 
used, we assume that the sampling rate of the signal is larger than any of the time delays due to the 
radar depth of the target. The magnitude squared of the Fourier transform of Ep in a frequency 
band (/, / + df) is proportional to the total radar cross-section <Tp due to many scattering centers 
located in a strip of width dY parallel to the X-axis centered at Y corresponding to f . The specific 
radar cross-section from this strip is, 

a° P {Y) = — f = / crp(X ,Y) sec 0 dX (11) 

d\ y_yfiTZyT 

where the sec 0 = R/y/R 2 — X 2 — Y^ factor is due to the fact that an area dA = R? sin © d© d$ 
located at polar angle © on the spherical surface projects to an area s ecQdX dY on the XY-plane. 
The total radar cross-section of the target is then given by one more integration, 

r R 

<t p = <T° P (Y)dY (12) 

J-R 

Once we calculate the average cross-section of a single backscattering event, that is, (<rp(X, Y)), 
then we can calculate the specific and total cross-sections, polarization ratios, and Doppler power 
spectra to compare with the data. In the next section we calculate this for the radar glory model. 


3.2 Backscattering from Buried Craters 

Figure 5 shows a buried crater from which we wish to calculate the radar cross-section. To analyze 
the backscattering we use geometrical optics to trace the rays through the ice overburden, through 
multiple bounces in the crater, and back through the ice-overburden. In the geometrical optics 
approximation the effects on the amplitude and field divide into purely geometrical considerations 
relating to the electrical path length of the rays and the spreading of a bundle of rays into a curved 
wavefront due to the curvature of the boundaries encountered and purely plane wave refraction and 
reflection effects at the boundaries (the usual assumption is that the radii of curvature of boundaries 
are large compared to the radius of curvature of the wavefront). Figure 5 shows the multiple effects 
on the complex amplitude of the electric field. First, the incident Electric field E x is multiplied 
by a rotation matrix J to put its components into TE and TM orientations. Then this field is 
multiplied by the Fresnel transmission coefficients for TE and TM incidence which is accomplished 
by multiplying by the matrix T tn . Then we multiply by a rotation matrix K to put this transmitted 
field’s components into TE and TM orientations for reflection from the walls of the crater. It 
turns out that for a spherical geometry such as we are assuming for the crater walls, the multiple 
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bounces lie in a plane, so the TE and TM components do not change orientation at each bounce. 
Then at each bounce we multiply by the Fresnel reflection coefficients appropriate for TE and TM 
incidence contained in the matrix R . These coefficients can be complex valued. Then we multiply 
by a geometrical divergence factor D accounting for the fact that after reflection from the curved 
boundary, i.e., the crater walls, the wave front’s curvature is altered. Since the geometry at each 
bounce is the same, after N bounces through the crater we multiply by N reflection matrices R N and 
the divergence factor D w (it is not simply a product of N divergence factors). Then, we multiply by 
a rotation matrix AT -1 to put the electric field components into TE and TM orientations appropriate 
for refraction back out of the ice which then contributes the factor T out • Then we multiply by the 
divergence factor D oui appropriate for transmission of a curved wavefront through a planar interface. 
Finally, we multiply by the rotation matrix J~ l to put the electric field components back into the 
orientations of the incident field. 

It turns out that there are a few other effects that must be considered. Figure 5 only shows one 
possible ray that makes it back to Earth. There are an infinite number of rays possible that make 
it back to Earth. If we look below the surface of the ice at the intersection of all these rays with a 
plane orthogonal to their path, we find they all lie on a circle of radius r 0 sin 0 where r 0 is the radius 
of the crater and 0 is the angle between the incident ray and the normal of the crater at each bounce 
in the crater. If we look in the plane orthogonal to the Earthward bound rays above the ice (shown 
in Figure 5), this circle becomes an ellipse in general. So, to compute the total field at Earth we 
have to add up the contribution from all these rays. The result is infinite; there is a caustic in the 
geometrical optics field in the backscattering direction. A number of methods are available to correct 
for this result without totally abandoning the geometrical optics approach. One method is to use the 
geometrical optics field on the plane indicated in Figure 5 and the Kirchoff diffraction integrals in the 
Fraunhoffer approximation to propagate the fields to Earth. Another method is Maslov's method 
(Maslov and Fedoriuk, 1981; Ziolkowski and Deschamps, 1984; Kravtsov, 1968) which constructs 
the geometrical optics field in 6 dimensional phase space (3 spatial and 3 wavevector components) 
where there is no caustics anywhere and then uses a hybrid 3 dimensional space consisting of some 
spatial and wavevector dimensions to compute the field on the caustic. Either method gives the 
same result which is finite. 

Another effect is that at some value of 0 the crater rim will obstruct some of the rays and prevent 
them from returning to Earth. The result is that the glory circles below the ice, ellipses above the 
ice will not be complete; they will break up into two arcs. Another effect is that we don’t expect 
that these craters are perfectly shaped and it is proposed that these glory arcs break up into patches 
or glints (Figure 6). The fields from each glint are assumed to be coherent, but the fields from two 
different glints add incoherently. The one exception is the two glints that are at opposite ends of 
a diameter of the glory circle below the ice. For a ray that travels through the path indicated in 
Figure 5, there is another ray that travels the exact opposite path. The fields associated with these 
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two rays add coherently and they come from the two glints mentioned. 

The radar cross-sections for fields received in the orthogonal and same linear (OL and SL respec- 
tively) polarizations and those for fields received in the orthogonal and same circular (OC and SC 
respectively) polarizations are 


<*«.(©, *) = <M t (3-* L ) , <r OL (e,*) = C7Ul + * L ) 
<7sc(0) = 2CM c (l + x c) , a oc {Q) = 2CA c (\-x c ). 


C- 


Ao 


(13) 

(14) 

(15) 


where F is the fraction of the glory arcs filled by glints, H is defined in Figure 6, and A 0 is the 
radar wavelength. The A's and x’s are (assuming total internal reflections occur; I’ll write the more 
general expressions later) given by 

A L = ^o{[(T? + T 2 ) + (T 2 - T 2 ) cos2($ - *p)p o - (1 + x)S] + 2T 1 T 2 (1 + x)5} (16) 


A l X l = ao{ 2 TiT 2 x<^o 

+[(T 2 - T 2 ) cos 2(<I> - <*> P ) - (T 2 + T 2 ) cos 4(<I» - * P )]fa> - (1 + *)tf] 
+2Ji7j[3(l + x)S - x<t> 0 ] cos 4(4* - <* P )} 

= ^o{(T 2 + T 2 )^o - (T\ - T 2 ) 2 ( 1 + x)6} 

A C X C = 2<toT 1 7Wo 


where 




sin 20 

V 

N 


7i o l sin 0' = sin©, is the direction of polarization for the linear polarization case, 


T\ = 4 oi cos 0 cos 0'/(cos 0 4- n 0 i cos 0') 2 


(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 


is the product of the TE Fresnel transmission coefficient for the ray entering the ice and the TE 
Fresnel transmission coefficient for the ray exiting the ice, 

T2 = 4tioi cos 0 cos 0' /(riQ\ cos 0 4- cos 0') 2 (22) 


is the same product for the TM fields, 

x = cos[27V tan -1 (cos 0(sin 2 0 — n 2 ) 1 ^ 2 / sin 2 0)] (23) 


is the mode decoupling factor defined in Eshleman (1986), 


6 = 0 0 /4 - sin(4^ 0 )/16 , 


(24) 
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and <j>o = 7r/2 if no shadowing occurs or sin</> 0 = cot 0 cot ©' when this last definition has a real 
solution. The factor D accounts for possible defocusing of the wavefront as it emerges from the ice 
into vacuum. Since Eshleman (1986) doesn’t consider the effects of the ice overburden he doesn’t 
include this factor, i.e., he has 

V = 1 . (25) 

If the ice overburden is very thick, then the wavefront encountering the ice is spherical and 


V = 


cos 2 0 
IXq X cos 2 0' 


(26) 


. We believe these two formulas for V provide upper and lower bounds respectively on the cross- 
sections. 

Figure 7 shows the dependence of A and x on the relative refractive index, n, of the layer of 
lower refractive index below the crater (relative to the refractive index of the ice overburden) and 
the number of bounces, N, for a crater at © = 0. The best fit to the observed total power in the 
various polarizations is given by x = 0.23 and N = 3. 

Figure 8 shows the variation in the radar cross-sections with respect to colatitude 0 and the 
shapes of the glory arcs. Surface craters and buried craters are compared. Of course the size of the 
radar cross-sections for the surface craters are much smaller than those of the buried craters. 

If we assume the density of craters is uniform across the surface of the satellite, then we can use 
eqs. (10), (11), and (12) to compute the Doppler spectra and total cross-sections. The resulting 
equations would look similar to those given above for a single crater except that now the A 7 s and 
x’s are integrals of the previous expressions across the surface and the constant C becomes, 


_ M 7rr a 2 FH 
4? iB? X 


(27) 


where M is the total number of buried craters on the surface of the satellite. Thus, C is proportional 
to the density of buried craters on the surface and is a free parameter to be determined by fitting 
to the observations. 

(Eshleman, 1986b) noted that a wavelength dependence that has been observed in the relative 
amounts of power reflected in the orthogonal polarizations could be due entirely to a wavelength 
dependence in the parameter x. One way that x may be wavelength dependent is if the refractive 
index at the crater walls makes a smooth transition from n 01 to n 02 . Here we will model the transition 
zone as a locally planar layer of thickness d with an exponentially decreasing refractive index (Figure 
9) because this is the simplest profile that admits an analytic solution for both TE and TM modes. 
The solution is in terms of Bessel functions whose order depend on kl sin 0 where / is the scale height 
of the exponential and 9 is the angle of incidence of the rays at the crater wall. 
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3.3 Comparison of Theory with Observations 

In the previous papers (Eshleman, 1986a; Eshleman, 1986b; Eshleman, 1987), it was found that 
x=0.23 is required to achieve a close fit to the observations for all three moons. There are several 
ways to get £=0.23, but the one that seems most likely is the three-bounce case involving TIR. 
This requires the hypothesis that n=0.7967. The first three rows of Table 6 compare the observed 
total radar cross-section and the ratios of the orthogonal components of the circular and linear 
polarization cross-sections with the predictions of the model. For & under the theory column we 
have used (7=1 and given the upper and lower bounds based on the new analysis and the analysis 
in Eshleman (1986) respectively. In order to match & for Europa we require (7=3.9 or 0.7 depending 
on the importance of a defocusing effect that was included in the above formulas. If FHf A were 
10 (say F=.5, H=2m, and A=.lm), then the buried crater model requires 39% or 7% coverage of 
the moon’s surface by glory craters. Also, note the excellent fit between the polarization ratios for 
Europa and Ganymede at 12.6cm. 

The last two rows of Table 6 characterize the shape of the power spectra in terms of an exponent 
for a cos m 0 fit to the actual curves. The power in narrow frequency bands is an integration over 
many craters on lines which are orthogonal to the line of sight and parallel to the spin axis of the 
moon. The theoretical numbers are based on the new analysis that we have called the lower bound. 
We don’t know yet to what degree the depth of the glints below the surface will affect the shape of 
the spectra. A comparison between the m values provide new support for the buried crater model. 
However, they do indicate that the theoretical spectra are narrower than the observed spectra. Also, 
the m oc /m sc values indicate that the observed OC spectra are a little narrower than the SC spectra 
while the model has the opposite trend. 

It is not clear how precisely a theory needs to match these exponents since the observed power 
spectra are full of additional structure that is not modelled well by the assumed law. The plots in 
Figure 10 show typical spectra for Europa, the cos m 0 fit to the data, and the spectra computed 
for the buried crater model. This clearly shows that a precise fit to m is not necessary to produce 
simulated spectra that model the data quite well. We also show for comparison simulated spectra 
for Earth’s moon based on empirical formulas determined by Evans (1957), where the SC echo is so 
weak that it is not visible at the scale shown. 

In Table 5 we see that Ganymede has been studied at two wavelengths and that there is a wave- 
length dependence in these observations. Here, we have only addressed the dependence indicated 
for <7 sc /(7oc • To explain the /i c data for Ganymede, x must go from about 0.22 at A=12.6 cm to 
about 0.33 at A=3.5 cm. Based on the assumption of a gradient in the refractive index at the crater 
walls as described in the last section, we have computed, in Figure 11, £ as a function of d/A for a 
single crater at 0 = 0 . We have found that the x value of a single crater is very nearly the same as 
the £ value for the entire moon. The 12.6 cm point is determined where £=0.22 and then the 3.5 
cm point is determined by the ratio of wavelengths. The x value is seen to increase with decreasing 
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wavelength which is consistent with the data. If this is indeed the cause of the dependence, radar 
measurements at different A would be extremely sensitive to the specific subsurface conditions on 
these moons. 



Chapter 4 


Conclusion 


The radar glory model is able to account for all the main features in the radar backscattering data 
from the icy Galilean satellites. In the previous chapter we presented the geophysical model in terms 
of a special structure involving buried hemispherical craters in which the refractive index abruptly 
decreases at the crater wall. It is not known if such a structure actually exists or how it could be 
formed in detail. It may even be possible that the structure is not formed from a cratering event 
and that the required variation in refractive index is due to changes in the ice crystalline state 
as discussed in Section 2.1. We were led to this model via the electromagnetic requirements for 
producing the observed echoes. It does not seem possible to produce the large radar cross-sections 
with a model such as that of Hapke (1990) or of Goldstein and Green (1980) involving complete 
randomness in the scattering mechanism. It seems that something more than the factor of two 
enhancement due to the coherent backscattering opposition effect discussed by Hapke is needed and 
that is what leads us to believe that a low degree of randomness imposed on a basically deterministic 
geometry that favors backscattering is required. In addition to this, though, total internal reflection 
also needs to occur to produce the large cross-sections and this conclusion is common to both the 
radar glory model and the Goldstein and Green model. Thus we conclude that these main features of 
the electromagnetic model — the low degree of randomness imposed on a deterministic geometry that 
favors backscattering such as a hemisphere and total internal reflection — are important ingredients 
of the true explanation and provide some information regarding the geophysical structures required. 
Unfortunately, the data do not sufficiently constrain the models and additional data types are needed 
to discriminate between the competing theories. 
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Table 1: 

Summary of the Physical Characteristics of the Icy Galilean Satellites 


Parameter Europa a Ganymede a Callisto a 


Orbital Radius b (10 3 km) 

Radius (km) 

Rotational Period (hours) 

Density (10 3 kg m“ 3 ) 

Subsolar Surface Temperature 0 (K) 
Average Surface Temperature 0 (K) 
Geometric Albedo 


670.9 

1,070 

1,883 

1,569 

2,631 

2,400 

85.22 

171.7 

400.5 

2.97 

1.94 

1.86 

132 

155 

169 

93 

110 

120 

0.6 

0.4 

0.2 


a Unless otherwise stated, data are from Burns (1986). 
b The orbits of all these satellites are virtually circular. 
0 Data from Gaffney and Matson (1980) 


Table 2: 

Absorption Length in Ice at Cold Temperatures a 
n = 1,8 


T (K) 


1/a (km) 


3.5 cm-A 

12.6 cm-A 

70 cm-A 

90 

7.1 

66. 

1.3 x 10 3 

no 

0.94 

8.8 

1.7 x 10 2 

130 

0.23 

2.2 

43. 

150 

0.08 

0.78 

15. 


Based on extrapolation formulas given by Thompson and Squyres 
(1990). 
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Table 3: 

Polymorphs of Ice That May be Present on the Icy Galilean Satellites 


Satellite 

Polymorphs stable 


at equator 


Europa 

II, VI, V, IX 

Ganymede 

II, Ic 

Callisto 

Ic 


Based on results from Gaffney and Matson (1980) 


Table 4: 

Bandwidths and Time Dispersions for the Icy Galilean Satellites 


Target 

T b (ms) 


B a (Hz) 



BT C 


3.5 cm- A 

12.6 cm-A 

70 cm-A 

3.5 cm-A 

12.6 cm-A 

70 cm-A 

Europa 

10 

3670 

1020 

184 

36.7 

10.2 

1.84 

Ganymede 

18 

3060 

850 

153 

55.1 

15.3 

2.75 

Callisto 

16 

1190 

330 

59.4 

19.0 

5.28 

0.95 


a Bandwidth, B = SirR/(\oP) with data from Table 1. 
b Time Dispersion, T = 2 R/c with data from Table 1. 
c A target for which BT > 1 is said to be overspread and it is not 
possible to simultaneously map the surface of the target in Doppler 
and time delay. 
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Table 5: 

Summary of Radar Data from the Icy Galilean Satellites 



Europa 
12.6 cm 

Ganymede 
12.6 cm 3.5 cm 

Callisto 
12.6 cm 

Earth’s Moon 
10. cm 

<7 

2.60 ± .64 

1.52 ±0.4 1.20 ±.28 

0.64 ±.16 

.05 

Pc 

1.56 ± .11 

1.55 ±.06 2.0 ±.03 

1.19 ±.06 

.05 

PL 

0.47 ±.07 

0.47 ± .08 

0.55 ±.10 

< 1 . 

m 

1.73 ± .08 

1.46 ± .04 

1.43 ±.05 

30. 

m 0 c /msc 

1.08 ±.10 

1.21 ± .06 

1.13 ±.09 

20. 


The data for the Galilean satellites are from Ostro (1982). Moon data are from Hagfors and Evans 
(1968) and Evans (1957): <7 is at 10 cm; the other numbers are at other wavelengths where they 
were available. 

& = (radarcross — section )/(ttR 2 ) 

<t = er sc 4- d’oc = <x SL + ^ol 
Me = vsc/voc 
Pl = vol/vsl 

m, m oc and m sc are exponents for a least squares fit of the power spectra to a cos m © law where 
0 is the angle between the radar line of sight and the surface normal of the moon. 

Where R is the radius of the particular satellite, C and L refer to circular and linear polarizations, 
and S and O refer to same and orthogonal to the transmitted polarization. 


Table 6: 

Comparison Between Observations and Buried Crater Model 



Europa 
12.6 cm 

Ganymede 
12.6 cm 3.5 cm 

Callisto 
12.6 cm 

Buried craters 
N=3, n=0.7967, C=1 

<7 

2.60±0.64 

1.52±0.40 

1.20±0.28 

0.64±0.16 

0.67 - 3.6 

^sc/d*oc 

1.56±0.11 

1.55±0.06 

2.0±0.03 

1.19±0.06 

1.60 

<7ol/ d'sL 

0.47±0.07 

0.47±0.08 


0.55±0.10 

0.44 

m 

1.73±0.08 

1.46±0.04 


1.43±0.05 

1.95 

m 0 c/ m sc 

1.08±0.10 

1.21±0.06 


1.13±0.09 

0.989 


The data for the Galilean moons are from Ostro in /7/. & = & sc 4- <r oc = <r SL 4- <r OL . m, moc and m sc are 
exponents for a least squares fit of the power spectra to a cos m 0 law where 0 is the angle between the radar 
line of sight and the surface normal of the moon. 
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Figure 1: AreciboGalilean satellites geometry. Geometry of radar backscattering from the Icy 
Galilean satellites. 


Figure 2: Satellite Geometry. Convenient coordinates for radar backscattering from a scattering 
center on a satellite. 


Figure 3: Comparison of radar echoes from a terrestrial planet and an icy Galilean satellite. 
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Figure 4: Five models for backscattering by icy satellites 


Figure 5: Buried crater geometry and its effects on radar backscattering 


Figure 6: Radar glory arcs and glints 
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Figure 7: Dependence of A and x on n and N 


Figure 8: Variation of buried and surface crater radar cross-section with colatitude 


Figure 9: Reflection of electromagnetic field from a layer with an exponentially varying refractive 
index 
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Figure 10: Comparison of buried crater Doppler spectra with observations from an icy Galilean 
satellite 


Figure 11: Variation of x for reflection from crater with a layer of material with exponentially varying 
refractive index 
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Fig 2. Typical power spectra for circularly polarized waves reflected from Europa compared 
with the theory and the c os" 1 0 law. T he data *re from O-stro 11 / page 22 1. The exponents 
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Developments in the Buried-Cra ter Model for Radar Echoes from lev Moons 
E.M. Gurrola, V.R. Eshleman (Stanford University) 

The buried-crater model introduced by Eshleman (Science 234 587, 1986) has 
been successful in accounting for anomalously large radar echoes with unex- 
pected polarization ratios that have been observed from three icy moons of 
Jupiter: Europa, Ganymede, and Callisto. A simple model of a single buried 
crater with axially incident rays plus a single assumption regarding the ratio 
of refractive indices in the layers above and below the crater walls was found 

adequate to model the power in two orthogonal circular and linear polarizations 
tor all three moons. z . q 

We have found that, with further development of 
the model to include effects due to angular po- 
sition of the craters plus a simple assumption of ' * 
uniform coverage by identical craters, this model 
also provides a good fit to the only other type 
of well-established data: the doppler-frequency 1 
power spectra. The figure shows theoretical curves 
of total radar cross-section for buried and surface 
craters and for a typical Europa spectrum (data •* 
from S. Ostro in Satellites of Jupiter, D. Morri- 
son, Ed., U. Arizona Press 1982). The light, un- 
labeled line is a least-squares fit of the data to a . # 
cos m © law where 0 is the angle between the 

incident rays and the moon’s surface normal. This fit to the data requires 8 times 
as many surface craters as buried craters. The buried-crater fit gives circular 
and linear polarization ratios (ratio of power in the unexpected to the expected 
polarization) of 1.60 and 0.44; averages of the observations for Europa are 1 56 
and 0.47 respectively. We also find by accounting for these additional effects, 
the original estimates of the reflected power from buried craters were too large.' 
However, the buried crater model still predicts larger echoes than other models 
that can also account for the polarization ratios. We will compare this model 
with other models that have been proposed with regard to their correspondence 
to current data as well as their predictions for other experiments that could be 
done in the future, such as range mapping. Despite an uncertain geophysical 
hypothesis, the buried-crater model is distinguished by its ability to simultane- 
ously account for the main features in all the types of data with a minimum of 
ad hoc assumptions. Finally, by introducing a gradation in refractive index at 
the crater walls, it is possible to model the sparse data indicating a frequency 
dependence in the polarization ratios. 
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ON THE ANGLE AND WAVELENGTH 
DEPENDENCIES OF THE RADAR 
BACKSCATTER FROM THE ICY GALILEAN 
MOONS OF JUPITER 

Eric M. Gurrola and Von R. Eshleman 

Center for Radar Astronomy, Stanford University, Stanford, CA 94305, 

U.S.A. 


ABSTRACT 

In this paper we report new developments in the buried crater model that has proved successful in ex- 
plaining t£ anomalous strengths and polarizations of the radar echoes from the rcy Gahlean moons o 
Jupiter - Europa, Ganymede, and Callisto. We have extended the theory to mate predictions of the radar 
cross-sections at all points on the surface of the moon, to compute the shape and strength of the power 
spectra, and to model a wavelength dependence that has been observed. 

DESCRIPTION OF THE MODEL AND ITS THEORY 

In 71-3/ Eshleman showed that the huge amplitudes and the unusual polarization properties of the radar 
echoes from Europa, Ganymede, and Callisto could be explained using a model of buried craters which 
have a decrease in refractive index at the crater walls. The electromagnetic model involves total intern'd 
reflection (TIR) in a geometry that strongly favors backscattering and is nnhr to the well known optuml 
phenomenon called the glory. This model has some features in common with the surface crater mode 
Ostro and Pettengill /4/, the random-facet model of Goldstein and Green /5/, and the refraction scattering 
model of Hagfors, Gold and Ierkic /6/. 

In the previous analysis, Eshleman considered the radar properties of a single buried crater at the center 
of the disc of each moon. This corresponds to 6=0 in the local crater geometry shown in Figure 1. Also, 
he did not enter into the complications due to the boundary at the surface of the ice. This interface has a 
small effect due to the plane wave transmission coefficients which may be obviated by a thin trumtion zone 
that matches the impedances of the two media above and below the surface. However this interface has a 
more significant effect in defocusing the curved wavefront that is returning to Earth from small cohe ™ 
glints in'the crater. This effect is dependent on the size of the two radn of curvature of the ^ fr ° 
relative to the depth of the glint. The previous analysis may be considered an upper bound on the modeh 
To bracket the results with a lower bound we have assumed that the glints are deep below the surface so 



Fig. 1. Geometry of a buried crater and a three-bounce (N=3) ray used in the electromag- 
netic model. The angle © is between the surface normal of the moon and the radar line of 
sight. For an ice-overburden, n«i = 1.8. In the model n«) < n*. (».e., n<l). 

(1)195 
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The simplest version of the theory obtains when the moon is assumed to be uniformly covered by identical 
craters. This may correspond in the actual situation to a uniformity in the average crater properties across 
the surface of the moon. In that case we can write the normalized (divided by *R r ) cross-sections in the 
following form : 

<7 sl = CA(3 — x) , 6oi = CA(l+x) 


£*c — 2CA(1 + x) , & oc = 2CA(1 — x) 

= Mirrj Fff 
4 xK* A 

where subscripts L and C refer to linear and circular polarizations; S and O refer to same and orthogonal to 
the transmitted polarization; M is the total number of buried craters on the surface of the moon; F is the 
fraction of the central glory circle that is filled by glints; H is the azimuthal size of each glint; r 0 and R are 
the crater radius and the moon radius respectively. The constant C is proportional to the surface density 
of craters and is used to model the observed radar cross-section of each moon. A and x are complicated 
integrals of the radar cross-sections of each crater over the surface of the moon. In the theory presented 
in / 1*3/, the mode decoupling factor, x, when TIR occured, was simply the cosine of the product of the 
number of bounces made by the rays in the crater and the difference in phase between the TE and TM 
reflection coefficients at the crater walls. In our new theory, x is more complicated in principle. However, 
numerically the x of the whole moon is virtually the same as the x for a single crater and virtually the 
same as in the original theory. 


The advantage of writing the cross-sections in the above form is that only two independent quantities - 
C.4 and x - are required to characterize four quantities - the four radar cross-sections of each moon. If 
the average crater properties are not sufficiently uniform then this is not possible; the x for the linear 
cross-sections will differ from the i for the circular cross-sections. Remarkably, the x values inferred 
separately for linear and circular polarizations indicate that Europa and Ganymede are consistent with 
the assumption of uniformity. 


In / 1 /, Eshleman noted that a wavelength dependence that has been observed in the relative amounts of 
power reflected in the orthogonal polarizations could be due entirely to a wavelength dependence in the 
parameter z. One way that z may be wavelength dependent is if the refractive index at the crater walls 
makes a smooth transition from n., to n„. Here we model the transition zone as a locally planar layer 
of thickness d with an exponentially decreasing refractive index because this is the simplest profile that 
admits an analytic solution. (The solution is in terms of Bessel functions whose order depend on kl sin 9 
where / is the scale height of the exponential and 9 is the angle of incidence of the rays at the crater wall.) 

RESULTS 


In the previous papers /1-3/, it was found that z=0.23 is required to achieve a close fit to the observations 
for all three moons. There are several ways to get x=0.23, but the one that seems most likely is the three- 
bounce case involving TIR. This requires the hypothesis that n=0.7967. The first three rows of Table 1 
compare the observed total radar cross-section and the ratios of the orthogonal components of the circular 
and linear polarization cross-sections with the. predictions of the model. For a under the theory column 
we have used C - 1 and given the upper and lower bounds based on the new analysis and the analysis in 
/ 1 / respectively. 



Europa 
12.6 cm 

Ganymede 
12.6 cm 3.5 cm 

Callisto 
12.6 cm 

Buried craters 
N=3, n=0.7967, C 

b 

2.60±0.64 

1.52±0.40 

1.20±0.28 

0.64±0.16 

0.67 - 3.6 

bsc/ Gqc 

1.56±0.11 

1.55±0.06 

2.0±0.03 

1.19±0.06 

1.60 

*Ol/ &SL 

0.47±0.07 

0.47±0.08 


0.55±0.10 

0.44 

m 

1.73±0.08 

1.46±0.04 


1.43±0.05 

1.95 

moc/nisc 

1.08±0.10 

1.21±0.06 


1.13±0.09 

0.989 


Table 1: The data for the Galilean moons represent averages over many observations and 
are from Ostro in jl / . b = <7 SC + ^oc = + b OL . m, m<>c and m sc are exponents for a 

least squares fit of the power spectra to a cos m 0 law where © is the angle between the 
radar line of sight and the surface normal of the moon. 


11 



Backscatter from Icy Galilean Moons of Jupiter 
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In order to match a for Europa we require C=3.9 or 0.7 depending on whether the lower or upper bound 
is used. If FH/A were 10 (say F=0.5, H=2m, and A=0.1m), then the buried crater model requires 39% or 
7% coverage of the moon’s surface by glory craters. Also, note the excellent fit between the polarization 
ratios for Europa and Ganymede at 12.6cm. 

The last two rows of Table 1 characterize the shape of the power spectra in terms of an exponent for a 
cos m © fit to the actual curves. The power in narrow frequency bands is an integration over many craters 
on lines which are orthogonal to the line of sight and parallel to the spin axis of the moon. The theoretical 
numbers are based on the new analysis that we have called the lower bound. We don’t know yet to 
what degree the depth of the glints below the surface will affect the shape of the spectra. A comparison 
between the m values provide new support for the buried crater model. However, they do indicate that 
the theoretical spectra are narrower than the observed spectra. Also, the m 0 c/m JC values indicate that 
the observed OC spectra are a little narrower than the SC spectra while the model has the opposite trend. 

It is not clear how precisely a theory needs to match these exponents since the observed power spectra 
are full of additional structure that is not modelled well by the assumed law. The plots in Figure 2 show- 
typical spectra for Europa, the cos" © fit to the data, and the spectra computed for the buried crater 
model. This clearly shows that a precise fit to m is not necessary to produce simulated spectra that model 
the data quite well. We also show for comparison simulated spectra for Earth’s moon based on empirical 
formulas determined by Evans /S/, where the SC echo is so weak that it is not visible at the scale shown. 



Fig 2. Typical power spectra for circularly polarized waves reflected from Europa compared 
with the theory and the cos m © law. The data are from Ostro /7 / page 221. The exponents 
for the data( theory) are: mo C =2.14(1.94), m sc =1.22(1.96), and m=1.56(1.95). The data 
and the theoretical curves have been normalized to <7=2.6; however the cos m © curve is 
not normalized in this manner. 


In Table 1 we see that Ganymede has been studied at two wavelengths and that there is a wavelength 
dependence in these observations. Here, we only consider the dependence indicated for £ S c/£oc- To explain 
the pL C data for Ganymede, x must go from about 0.22 at A=12.6 cm to about 0.33 at A=3.5 cm. Based 
on the assumption of a gradient in the refractive index at the crater walls as described in the last section, 
we have computed, in Figure 3, x as a function of d/A for a single crater at © = 0 . We have found that 
the x value of a single crater is very nearly the same as the x value for the entire moon. The 12.6 cm point 
is determined where x=0.22 and then the 3.5 cm point is determined by the ratio of wavelengths. The x 
value is seen to increase with decreasing wavelength which is consistent with the data. If this is indeed the 
cause of the dependence, radar measurements at different A would be extremely sensitive to the specific 
subsurface conditions on these moons. 
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0.00 0.05 0.10 0.15 0.20 0.25 

d/A 

Fig 3. x versus d/A for three values of n near the optimal value determined for A=12.6cm. 
CONCLUSION 

We have demonstrated that the buried crater model is consistent with the anomalous radar c h a r acteristics 
of the Galilean moons. It can explain the surprisingly large cross-sections, the odd reversal of expected 
polarizations, the broad angular spectrum of the returns, and the apparent change in polarization properties 
with radar wavelength. More accurate and more discriminating types of radar measurements are needed 
as further guides for determining the proper theory of the radar echoes from these moons. 
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Radar Backseat tering by The Icy Galilean Satellites and The Radar Glory Effect. 
Eric M. Gurrola and Von R. Eshleman. Center for Radar Astronomy, Stanford University. 


Radar echoes from the icy Galilean satellites — Europa, Ganymede, and Cailisto — are unique 
in character, and they may be explained by a phenomenon that has been called the radar glory 
effect. 1,2,3,4 We have analyzed the radar glory backseat tering from buried craters as a possible 
model for this effect. These craters have a smaller refractive index below the crater than above 
as shown in Fig. 1(A). The possibility exists that the rays shown will be totally internally 
reflected at the crater walls, and we assume this to be true in the following. The rays which 
contribute to the backscattering will come from a circular annulus, when viewed before the 
refraction occurs at the ice-vacuum surface, in a plane orthogonal to the rays reflected from 
the crater (see Fig. 1(B)). For craters near the limb of the satellite some of the rays may be 
blocked by the crater rim, and this circular annulus will break up into two circular arcs. We 
assume that these glory arcs are broken up into glints of dimension JT, assumed larger than the 
wavelength, and that the fields from a pair of glints connected by the ray shown add coherently 
while the fields from the various glint pairs add incoherently. The coherent addition of the fields 
from glint pairs arises because for the ray shown in Fig. 1(A), there is another ray that travels 
the exact same path, but in the opposite direction. 



Figure 1. (A) Geometry of a buried crater and a three-bounce (N=3) ray. The angle 
0 is between the surface normal of the satellite and the radar line of sight. For an 
ice-overburden, noi = 1.8. In the model n«a < noi («'•«., n<l). (B) Glory arcs of angular 
extent 2d>o and a pair of glints of dimension H. TE and TM electric field vectors are 
also shown. 


Consider a crater located near the surface of the satellite at coordinates (R,0,$), where R 
is the satellite radius, 0 is longitude with 0 = 0 in the direction of Earth, and * is latitude with 
$ = jt/ 2 in the plane defined by the spin axis of the satellite and the 0 = 0 direction. The radar 
cross-sections for fields received in the orthogonal and same linear (OL and SL respectively) 
polarizations and those for fields received in the orthogonal and same circular (OC and SC 
respectively) polarizations are 


*sL(e,*) = CA L (3-jr L ) , <r 0 L(e,*) = CA L (i + x L ) 
<r S c(0) = 2CAc(l + X c ) > 0bc(0) = 2CAc(l - A c ). 


( 1 ) 

( 2 ) 


The constant C is proportional to the number of glint pairs in the glory arcs, and is given by 
C = irrlFH/X where F is the fraction of the glory arcs filled by glints, H is defined in Fig. 
1(B), and A is the radar wavelength. The A’s and X’s are given by 

Al = \doWl + Tl) + (Ti - Tj )cos2($ - * P )]fa> - (1 + x)S] + 27\T 2 (1 + x)S) (3) 

= <7o{ 2Tj 

+[(I? - J?)cos2(# - ip) - (7? + 7?)cos4($ - $ P )][<&> - (1 + x)$] (4) 

+2Ti7a[3(l + x)6 — x0o]cos4($ — $ P )} 

A c = cr 0 {(Ti + Ti)<f > o - - T 2 f( 1 + x)6} (5) 

Ac Xq — 2<7 0 T] J 2 x (6) 

where 

. sin 20 cos 2 0 

°° N n^jcos 2 ©'’ ^ 

noi sin 0' = sin@, ip is the direction of polarization for the linear polarization case, T\ = 
4noi cos 0 cos 0'/(cos 0 + noi cos 0') 2 is the product of the TE Fresnel transmission coefficient 
for the ray entering the ice and the TE Fresnel transmission coefficient the ray exiting the 
ice, T 2 ss 4noi cos©cos07(noi cos© + cos©') 2 is the same product for the TM fields, x = 
cos [2^ tan -1 (cos0(sin 2 0 — n 2 ) 1|/2 / sin 2 0)] is the mode decoupling factor defined in reference 1, 
6 = <fo/4 — sin(4$o)/16, and fo = x/2 if no shadowing occurs or sin ^0 = cot 9 cot 0' when this 
last definition has a real solution. 

The power spectra are sums of the contributions from several craters along lines of constant 
Doppler shift which are lines, when projected onto the plane of the sky, orthogonal to the line 
of sight and parallel to the spin axis. That is, the cross-sections at a specified frequency are 
integrals of the cross-sections in equations (1) and (2) on semi-circles on the satellite surface 
where cos 0 cos i =constant. In Fig. 2 we compare the computed cross-section for circular 
polarization with echoes obtained from Eurpoa; this spectrum is typical of the spectra obtained. 

The total cross-section of each satellite is an integration of the cross-sections in equations ( 1) 
and (2) over the entire illuminated surface. These total cross-sections can be written in a form 
similar to those in equations (1) and (2) except that the j4’s and X’s are the same for linear 
and circular polarizations and are complicated integrals that must be computed numerically. 
In Table 1 we compare a summary of the main features of the data with the model predictions. 
The scaling constant C for these total cross-sections is given by C = ( Mrro/4irR 2 )(FH /\ ) 
where M is the total number of craters on the surface. 
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Figure 2. OC, SC, and OC+SC cross-sections obtained from Europa (dashed curves) 
on January 10, 1990 at the Arecibo Observatory at 12.6 cm-A and from the model (solid 
curves) using N = 3, n=0.7967. The data and the model have been scaled to give a 
total cross-section of 2.6. 


Europa Ganymede Callisto 


a 2.60 1.52 0.64 

dsc/doc 1-56 1.55 1.19 

d 0L /d S L 0-47 0.47 0.55 

m 1.73 1.46 1.43 


Buried craters 
N = 3, n=0.7967, C = 1 



0.44 

1.95 


Table 1: The data for the Galilean satellites are for 12.6 cm-A and are from Ostro 5 . 
a = <j sc + < 7 oc = d SL -I- <7 0 l . These cross-sections are normalized by the geometrical 
area *R 2 , i.e., a = ajvS?. Here m is the exponent for a least squares fit of the power 
spectra to a € 08 “ 0 law. 
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